Does a single zealot affect an infinite group of voters ? 
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A method for studying exact properties of a class of inhomogeneous stochastic many-body systems 
is developed and presented in the framework of a voter model perturbed by the presence of a "zealot" , 
an individual allowed to favour an opinion. We compute exactly the magnetization of this model 
and find that in one (Id) and two dimensions (2d) it evolves, algebraically (~ t~^^^) in Id and much 
slower (~ l/\nt) in 2d, towards the unanimity state chosen by the zealot. In higher dimensions 
the stationary magnetization is no longer uniform: the zealot cannot influence all the individuals. 
Implications to other physical problems are also pointed out. 
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Researchers have devoted much attention to the field 
of non-equilibrium many-body stochastic processes [ij. 
In particular the study of exact solutions of prototypi- 
cal models such as the voter model 0] has proven to be 
fruitful for understanding a large class of non-equilibrium 
phenomena 0. However, more realistic systems taking 
into account the important dynamical effects of inhomo- 
geneities, constraints and disorder (see e.g. 0, Q and 
references therein) are less understood. To gain com- 
prehension of these situations, exact results for systems 
modeling inhomogeneous situations in any dimension are 
desirable, but rather scarce y. 

In this work, motivated by these considerations, and 
with the voter model as a paradigm, we present tech- 
niques for computing exact properties (in any dimension) 
of a class of stochastic many-body systems with inhomo- 
geneities. 

The voter model is an Ising-like model where an "in- 
dividual" (or spin) associated to a lattice site r can have 
two different opinions (Xr = ±1 0- The dynamics of the 
system is implemented by randomly chosing one individ- 
ual and assigning to it the value of the spin of one of its 
randomly chosen nearest neigbors. In the voter model, 
the global magnetization is conserved and the dynamics 
is Z2 symmetric (invariance under the global inversion 
(Tr — CTr)- The importance of the voter model stems 
from the fact that it is one of a very few stochastic many- 
body systems that are solvable in any dimension and is 
useful for describing the kinetics of catalytic reactions 
0, IH 111 J in studying coarsening phenomena 0, Q and as 
a prototype model of opinion dynamics Q . 

For the sake of concretness, and without loss of general- 
ity, we specifically present our method and techniques in 
the framework of an inhomogeneous voter model where, 
to mimic in a simple manner the fact that a group of 
agents may have heterogeneous interactions, the conven- 
tional voter dynamics is supplemented by the presence 
of a "zealot" : a biased individual who favours one opin- 
ion. We study the effect of this perturbation by com- 
puting the exact long-time magnetization, which is no 
longer conserved, in dimensions d = 1, 2 and 3. In low- 



dimensions the zealot drastically affects the dynamics: 
the system evolves towards unanimity with the latter. 
The approach to the stationary state is algebraic in Id 
and logarithmically slow in 2d. In d > 3, the effect of 
the zealot is less dramatic and the local stationary mag- 
netization is a non-trivial function of the distance to this 
biased individual. These findings, although formulated 
in an "opinion dynamics" language, are, as pointed out 
hereafter, relevant to other physical problems. 

The inhomogeneous voter model that we study is de- 
fined on a hypercubic lattice of size (2L -I- 1)'^, where 
individuals, labelled by a vector r having components 
—L < ri < L (with i — 1, . . . ,d), may interact accord- 
ing to the usual voter dynamics. In addition, we now 
consider that a zealot, at site "0", tends to favour the 
diffusion of the -1-1 opinion via the interaction with his 
neighbors: the zealot is the only individual in the system 
allowed to change his state from —1 to -1-1 (with rate 
7 > 0) without regards to his neighbors, with whom he 
nevertheless interacts. According to the spin formulation 
of the model, the state of the system is described by the 
collection of all spins: S = {<Jr}. In this language, the 
dynamics of the model is governed by the usual voter 
model transition-rate 0, S S supplemented by a local 
term involving the local zealot's reaction. The spin-flip 
rate, Wr{S) = w^a^ — > — Cr), therefore reads 

Wr{S) - ^ (^1 - ^^r E ^r'^ + J (1 - '^O) <5r,0. (1) 

Here the sum on right-hand side (r.h.s.) runs over the 
2d nearest neighbors r' of site r and r = l//3d > de- 
fines the time scale. The probability distribution P{S, t) 
satisfies the master equation: 

j^P{S, t)^J2 t) - w,{S)P{S, t)] , (2) 

r 

where the state 5"" differs from S only by the spin-flip of 
CTr- With the master equation in the limit L ^ oo, 
the equation of motion of the local magnetization Sr{t) = 



2 



^^(TrP(S', reads: 



ArSrit)- jiSo{t)-l)dr,0. (3) 



Here Ar denotes the discrete Laplace operator: 

The last term on the r.h.s. of Eq.© is due to the effect 
of the zealot and only appears in the case where r = 0. 

An important consequence of Eq.Q is that the only 
possible uniform final magnetization is S'r(oo) — 1, i.e. 
the state favoured by the zealot. However, in d > 3, 
as shown hereafter, the stationary magnetization profile 
turns out to be non-uniform (but isotropic). 

Using the properties of the modified Besscl functions 
of the first kind, Ir{t) we obtain the formal solution 
of Eq.Q: 



Sr{t) = ^5k(0)n [e"'^*4.-r.(2/3i)] 



1=1 



+ 7 / dt'[l^Soit-t')]ri\e-^''''lr,i2pt') 
'0 i=i 



(4) 



To obtain an explicit expression for the magnetization we 
solve the self-consistent integral equation @ for r = 
and then plug the result back into For this purpose 
it is useful to denote the Laplace transform of a product 
of Bessel functions (multiplied by an exponential term) 



/r(s,/3) 



dt e-^* [e-^'"''lr, (2/3t) (2/3i)] . (5) 



This quantity can be rewritten in terms of Watson inte- 
grals, or "lattice Green-functions": 



d'^q 



-2q.r 



C^ttY s + 2P[d- cos q,Y 



(6) 



where q = (gi, . . . , qd) is a d— dimensional vector. 

Laplace-transforming Equation Q , and using the con- 
volution theorem, we obtain the following expression for 
the Laplace transform of the local magnetization: 



Sris) 



dtt 



7^r(s,/3) 



l + 7/o(s,/3) 



(7) 



where, for technical simplicity, we consider the zero initial 
magnetization state: Su{0) — 0. 

The exact expression for the long-time magnetization 
is obtained by Laplace-inverting the s —^ expansion of 
Eq.(|Zj) and paying due attention to the situations where 
the integrals ((HJ are divergent. In the sequel the norm of 
the vector r is denoted r. 

We first consider the one-dimensional problem. In this 
case, the quantity © reads f§|: 



Us,f3) = ir{s,f3) = 



.(8) 



It should be noticed that the s ^ behavior of ^ 
diverges as s~-^/^. 

Laplace- inverting the s — > expansion of Q , together 
with the expression of © for r = 0, we get (/?< oo): 



(9) 



In the more general case, where r > 0, the local mag- 
netization is obtained similarly from j?)) and ^ [in the 
limit s 0, with r-^s kept fixed]: 



Sr{t) ~ erfc 



2y^ 



(10) 



where erfc(z) is the usual complementary error function 
0. The expression (|10|l is valid for /3t 3> 1, and is par- 
ticularly useful in the scaling limit where both r ^ oo 
and t oo, but where the ratio r/^/j3i is finite. 

When r is finite (i.e. < r < oo ) and pt ^ oo, we 
obtain the following long-time behavior: 



Sr{t) = 1 - 



r + (2/3/7) 



(11) 



In two dimensions, the integral JBJ is also divergent in 
the long-time regime s —^ and therefore its main con- 
tribution arises from q"^ = + q2 ^ 0. In this sense, we 
first expand © for small s in the case where r = 0: 



/o(s,/?)— In 
s^o 47rp 



(12) 



which implies that the long-time behavior of the average- 
opinion of the zealot is 



Soit) - S'o(oo) 



1 



V 7 y 



(13) 



where 50(00) = 1. 

For the other individuals we proceed similarly and from 
©, with r 1, we have: 



(14) 



where Kq{x) is the usual modified Bessel function of the 
third kind (9|. Using the small argument expansion of 
such a Bessel function 0, we find that the long-time 
behavior in the region where (3t ^ r^ ^ 1 is 



Sr(t) — S'r(oo) 



In r^ 



(15) 



where the stationary magnetization corresponds again to 
the unanimous opinion, as in Id: S'r(oo) = 1. 



The regime where r oc /3t 



should still be dis- 



cussed separately. In fact, in this regime, from IT)), l|12() 
and H14|l we have : 



2Kn 



Sr{s) 



sin {s/ 13) 



(16) 
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where s and r — > oo, but r^J s/ (3 is finite. 

From results (0, lEl) and lUll), CH) wc infer that 
in low-dimensions, at large time, the effect of the zealot 
appears at two length scales: (i) The opinion of individ- 
uals "close" (r^ <C fit) to the zealot evolves algebraically 
(oc t~^/'^) in Id and logarithmically (cx 1/lnt) in 2d to- 
wards the unanimous opinion 5r(oo) — 1. (ii) For indi- 
viduals "far" (r^ (x (3t ^ oo) from the zealot, the local 
magnetization evolves as a smooth scaling function of 
u = r^/2v^ in Id. This is however no longer the case 
in 2d, where, due to logarithmic terms arising in (|15|l . 
the magnetization has not a scaling form. A qualita- 
tively similar result, but within a different context, has 
recently been reported in Reference ^o|- can also 

notice that results , and (^HJ , IIEI) are (mainly) 
independent of the rate 7: for the long-time behavior of 
the local magnetization; in low-dimensions, only the fact 
that there is a biased individual (i.e. 7 > 0) matters. 

We therefore conclude that both in Id and 2d the zealot 
eventually affects all the individuals: the number of vot- 
ers having a final -1-1 opinion are within a "circle" whose 
radius increases as 00. 

We now study the three-dimensional situation, and 
then consider its generalization to the case where c? > 3. 

When d > 3, the integrals © are well defined for all 
the values of s, and in particular when s — )■ 0. There- 
fore, conversely to what happens in Id and 2d, to de- 
termine the long-time behavior of the magnetization we 
cannot simply focus on the (7 — > expansion of For- 
tunately, very recently Glasser and Boersma have been 
able to explicitly compute © in the 3d case where s = 
[Tl| . We now take advantage of these findings to compute 
the stationary magnetization in 3d. We therefore intro- 
duce a triplet (ar,6r,Cr) of rational numbers depending 
on r, given in Table 2 of Reference and the quan- 

tity 50 = (^) r2 {^) r2 (li), where T{z) is Euler's 

Gamma function ||9j]. 

With help of the results obtained in , we are in po- 
sition to compute the 3d stationary local magnetization 
(SLM) by taking the s ^ limit in 0: 



S'r(oo) 



go [2(3 + jgo] 



(17) 



In particular, for r = we have the triplet (oq, 60, co) = 
(1,0,0) 113 and thus obtain 5'o(oo) = 2^^- 

To gain, in a simple manner, some further insight of the 
behavior of the discrete expression (|17(l when r ^ 00, it 
turns out to be fruitful to take the continuum limit of . 
In this limit the SLM is 5(r,oo), and we have to solve 
the problem of determining the electric potential due to 
a "charge" at the origin. One should pay due attention 
to the fact that this electrostatic reformulation needs to 
be supplemented by additional information since, from 
the continuum limit of Q , the charge is a priori an un- 
known quantity. To overcome this difficulty, in 3d, we 
take advantage of our knowledge of the discrete version 



of the problem and, with l|17|l . we compute the charge 
at site r = as ^[1 — 5o(oo)] = 2/3+730 ' therefore 
obtain the following 3d continuum stationary equation: 



A5(r,oo) 



27 



2/3 + 7.90 



(18) 



where A is the 3d-Laplacian operator and d^{r) denotes 
the 3d-Dirac delta function. 

The solution of Eq. (|18|) depends on 7 and reads: 



5(r, 00) 



7 



1 



2^(2/3 + 730) r 



(r > 0). (19) 



Comparing the predictions of the results (fTTIl and ifT^ . 
we notice that they agree very well, even for finite r: For 
instance, when 7 = 2/3, at r = (3, 1, 0), we have 5'r(oo) = 
0.0339 and 5(r,oo) = 0.0334, whereas for r = (3,1,1) 
5r(oo) = 0.0319 and 5(r,oo) = 0.0319. 

Results (|17|l and 119|l show that in 3d (conversely to 
what happens in low-dimensions) the SLM is an isotropic 
(which is clear from ((HJ and (|19l) ) but non-uniform func- 
tion decaying with the norm of r: Sr{oo) « 5(r,oo) = 
^-i^ll^ where the amplitude is given by (|19|l [and r > 0]. 
The reasoning can be extended to dimensions d > 3, 
where the electrostatic reformulation gives the result 
(for r > 0): 5(r,oo) = ^r^ry. Again the computation 
of the amplitude Ad requires the explicit knowledge of 
/o(s = 0, /3) in dimensions d > 3. 

Despite the fact that in d > 3 the s ^ analysis of © 
is a difficult task, we can infer (for r finite) the long-time 
behavior from Sr{t) - 5r(oo) (/3i)^"^. 

From Eq.|01) we can also compute the total magneti- 
zation: M{t) = J2r Sr{t) = M(0) + 7 /o dr [1 - So{t)]. 
This expression shows that in the voter model (7 = 0) 
the quantity M{t) is conserved, which is no longer the 
case when 7 > 0. With the help of the results ®, l(T^ 
and (|17|l . we obtain the long-time behavior of AI (t) [here 
M(0) = 0]: 



Mit) 



(/3t)V2 
Pt/ Int 
(3t 



d = 1 
d = 2 
d > 3. 



{pt 00) (20) 



These results show that in this inhomogeneous voter 
model the saturation time t^, that is, the time necessary 
to have M{ts) comparable to the size L'^ (where L — > 00) 
of the system, scales as ptg ~ in Id, ptg ~ L^lnL 
in 2d, and pts ^ L"^ tor d > 3. These statements are 
in agreement with results obtained for other models (see 
e.g. 4]).^ 

In this work we have developed a method to compute 
some exact properties of a class of stochastic many-body 
problems with inhomogeneities and have explicitly pre- 
sented this approach in the framework of an inhomoge- 
neous voter model where the usual voter dynamics is per- 
turbed by the local presence of a single zealot. For this 



4 



"opinion dynamics" problem, we have computed exactly, 
in dimensions d — 1,2 and 3, the long-time magnetiza- 
tion (mean- average opinion of each voters). From our 
exact results we have seen that in low-dimensions the 
zealot (i.e. the inhomogeneity) always affects the mean- 
average opinion and that its effect propagates as t^^^. In 
fact the mean opinion of individuals approaches, alge- 
braically in Id [see (jllf) ]. according to the scaling expres- 
sion H10|l . and logarithmically slowly in 2d [see (|15|) ]. the 
unanimity state favoured by the zealot. These results are 
(mainly) independent of the strength of the biased indi- 
vidual. In 3d the situation is completely different and 
the (stable) stationary mean-average opinion of a voter 
is no longer uniform but follows the non-trivial isotropic 
function H17(l that decays with the inverse of the distance 
to the zealot [see (|19|l ]. a result also obtained in the con- 
tinuum limit via a suitable electrostatic reformulation, 
and then extended to the case where d > 3. The findings 
obtained in this work, and the differences between the be- 
haviors observed in low-dimensions and in d > 3, can be 
qualitatively understood in realizing that the local mag- 
netization is the solution of a diffusion-like equation ^ 
supplemented by a local boundary term (the zealot) and 
taking into account the fact that in Id and 2d random- 
walks are recurrent (which, in our case, implies that all 
individuals interact with the zealot), while in d > 3 they 
are transient, and therefore there is a finite probability 
that individuals never interact with the zealot Q, Q . It 
is instructive to compare our results to those obtained in 
the conventional voter model .Si, 13 • the presence 
of the zealot clearly implies that the magnetization is no 
longer conserved and that the dynamics is not transla- 
tionally invariant. This comparison also shows that a 
single inhomogeneity (here, the zealot) can deeply affect 
the stationary and the long-time properties of an inter- 
acting spin system, whose perturbed dynamical behavior 
depends on the dimension d. 



Despite the fact that our method and results have been 
formulated in an "opinion dynamics" language, we em- 
phasize that they have a broad physical relevance and 
can be applied to a large class of stochastic many-body 
problems. As a physical illustration we can consider the 
kinetics of the monomer-monomer catalysis on an inho- 
mogeneous substrate that locally dcsorps preferentially 
one species of monomer 4, 5, 12J. Identifying spins -1-1 
(— 1) with A [B) particles, the monomer-monomer re- 
action can be mapped onto an Ising model with mixed 
voter and Kawasaki dynamics 0, 01, whose time-scales 
are respectively defined by l/(3d and l/(3'd, and an in- 
homogeneous term [as in Q] that mimics the local des- 
orption (at r = 0), with rate 7, of particles of species B 
H [12. For this model, with /3 = /3 -f the local con- 
centration Cr{t) of A particles obeys to {d/{Pdt))cr{t) = 
ArCr(t) — (7//3)(co(t) — l)5r,o, whcrcas the concentration 
of B particles, at site r, is 1 — Cr{t) Comparing the 
equation for Cr with Eq. Q, it is clear that the concen- 
tration of particles in the catalysis problem can imme- 
diately be inferred from the above results for Sr{t) [l^. 
Another field of pertinence of this work, according to the 
well-known relation between the (Id) voter and Glauber- 
Ising models , is the area of inhomogeneous mag- 

netism. We can also mention that our method, which can 
take into account the presence of many inhomogeneities 
and can deal with systems of more than two states per 
site, is of direct relevance for a large class of reaction- 
diffusion models 0|. In this context, a diffusive model of 
interacting particles where an homogeneous source is in 
competition with a local trap has been solved 12]. 
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